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In the universality class of the one dimensional Kardar-Parisi-Zhang surface growth, Derrida and 
Lebowitz conjectured the universality of not only the scaling exponents, but of an entire scaling 
function. Since Derrida and Lebowitz's original publication [PRL 80 209 (1998)] this universality 
has been verified for a variety of continuous time, periodic boundary systems in the KPZ universality 
class. Here, we present a numerical method for directly examining the entire particle flux of the 
asymmetric exclusion process (ASEP), thus providing an alternative to more difficult cumulant ratios 
studies. Using this method, we find that the Derrida-Lebowitz scaling function (DLSF) properly 
characterizes the large system size limit (N — * oo) of a single particle discrete time system, even 
in the case of very small system sizes (N < 22). This fact allows us to not only verify that the 
DLSF properly characterizes multiple particle discrete-time asymmetric exclusion processes, but 
also provides a way to numerically solve for quantities of interest, such as the particle hopping flux. 
This method can thus serve to further increase the ease and accessibility of studies involving even 
more challenging dynamics, such as the open boundary ASEP. 

PACS numbers: 02.50.-r,05.40.-a,82.20.-w,89.75.Da 



X 



I. INTRODUCTION 

The Kardar-Parisi-Zhang (KPZ) equation [1] describes 
a rich variety of processes such as surface growth [2, 3], 
directed polymers [4-6], and avalanches [7, 8]. And ac- 
cordingly, the massive array of literature (see [9-11] and 
references therein), on the topic of KPZ theory reflects 
the central role the KPZ equation plays in the study of 
stochastic dynamic processes. The wide variety of dy- 
namic processes governed by the KPZ equation form the 
so-called KPZ universality class — a class of seemingly 
unrelated dynamics whose bulk properties obey, on a 
course-grained level, this one master equation. One mem- 
ber of the KPZ universality class, the well studied one- 
dimensional asymmetric exclusion process (ASEP) de- 
scribes a driven lattice gas with hard core exclusions [12], 
and has also been applied to studies of highway traf- 
fic [13, 14], protein synthesis [15], and sequence align- 
ment [16-18]. Characterizing the underlying properties 
of the ASEP promises a greater understanding for these 
specific studies as well as insights into the broader KPZ 
universality class. Not surprisingly, much effort has been 
spent calculating some of the many properties of the 
ASEP such as the density profile, steady states, mass 
gaps, and diffusion constants [19-28]. Though many 
questions still remain, these myriad studies have helped 
to uncover a great deal of insight into the ASEP. 

In their study of the ASEP, Derrida and Lebowitz [29] 
extended the Bethe Ansatz approach of Gwa and 
Spohn [23] in order to solve the totally ASEP for parti- 
cle displacement in the asymptotic limit of large system 
size [29] . One of the most interesting aspects of their SO- 
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lution involves the scaling function G, describing the non- 
linear behavior of the total particle flux. The Derrida- 
Lebowitz scaling function (DLSF) is independent of any 
model parameters. Thus, this scaling function was con- 
jectured to be universal, i.e., characteristic of all KPZ 
systems. In a follow-up study, Derrida and Appert an- 
alytically continued the Derrida-Lebowitz scaling func- 
tion (DLSF), successfully completing the solution for all 
space in the asymptotic limit of large system sizes [30]. 
Since then, a large number of studies have given strong 
evidence for the universality of the DLSF. These stud- 
ies fall into two classes. On the one hand, for a few 
closely related variants of the continuous time totally 
ASEP [4, 7, 31] the characteristic DLSF behavior has 
been analytically verified. On the other hand, numeri- 
cal studies have bolstered the universality claim of the 
DLSF [7, 30-33] for a much broader range of systems in- 
cluding genuinely discrete time systems. However, these 
numerical methods do not directly verify the universality 
of the DLSF, but rather verify the universality of certain 
cumulant ratios which must be universal if the DLSF is 
universal [30]. These numerical approaches cannot di- 
rectly verify the DLSF since they use sampling methods 
that are inherently unable to probe the full DLSF which 
contains information about statistically rare events. In 
addition to that, both analytical and numerical studies 
rely on examining the behavior of the so-called scaling 
region [30] which applies only in the limit of very large 
system sizes (~ 1280), especially in the discrete time case 
(> 10240). 

By examining properties applicable to the intermediate 
scaling region, we created a method for directly measur- 
ing the DLSF for discrete time systems of considerably 
smaller size (< 22). Our method works without resorting 
to the stochastic sampling that makes cumulant meth- 
ods so time consuming. As an application, we show that 
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FIG. 1: Diagram of a periodically bound totally ASEP. Only 
the particle shaded in gray, on site 3, may move to the space 
on the right, outlined by a dotted line, with a probability 
dt. The remaining particles are prevented from moving by 
an adjacent particle occupying the space to the immediate 
right. More explicitly, the particle on site 2 finds it's advance 
blocked by the particle that resides on site 3. Meanwhile, the 
particle on site 1 is likewise blocked by the particle on site 2. 
The particle on site N also finds it's movement thwarted by 
the particle on site 1 due to the periodic boundary conditions. 



the discrete time ASEP with both single and multiple 
particles per site under any parameter choice is charac- 
terized by the same DLSF conjectured to be universal 
for all processes within the KPZ universality class. Our 
method can also be used to study the particle hopping 
in other ASEP scenarios, including the open boundary 
ASEP. Once DLSF behavior has been verified for any 
given system, our method also provides a way to numer- 
ically calculate the non- universal scaling constants that, 
together with the known form of the DLSF, can be used 
to calculate properties such as the particle hopping dis- 
tribution or the large deviation function. 

This paper will begin with a short review of the con- 
tinuous time single particle system from which Derrida 
et al. [29, 30] initially derived the DLSF in section II. 
Then, we will reveal a small extension of the Derrida et 
al. results that will allow us to measure the DLSF for the 
discrete time ASEP. Next, we describe our method and 
its use with the discrete-time single- and multi-particle 
ASEP as an example of the application of our method. 
Last, we will discuss the applicability of this method to 
other systems of interest. 



II. REVIEW OF THE ASYMMETRIC 
EXCLUSION PROCESS 

Derrida et al. [29, 30] examine a periodically bound 
continuous time single particle totally ASEP with N sites 
each capable of holding one particle as shown in Fig. 1. 
These particles may move only to the right with a prob- 
ability dt if the target site is unoccupied. The total num- 
ber of particles in the system is fixed at p. 

In order to characterize the probability distribution 
P{yt) of the cumulative particle hopping per site yt, they 
calculate the generating function Z, given by 



Z = (exp[7y t 



(1) 



where the brackets (•) denote the ensemble average. In 
order to solve for this generating function, they use the 




FIG. 2: Plot of the left hand side of Eq. (3) rescaled by N 3/2 
versus yN 1 ' 2 for variously sized systems. A plot of the ap- 
propriately rescaled DLSF G is included for reference. The 
curves for smaller N do not agree well with the DLSF. This 
tells us that the left hand side of Eq. (3) cannot be described 
by G alone for small system sizes. The gap in the curves is a 
consequence of the dramatic increase in computational effort 
in its vicinity due to the slow convergence of the infinite series 
defining the solutions Aw (7) close to their radii of convergence 
(see appendix A). 



large time t behavior 



(exp[7j/ t ]) ~ exp[AAr(7)£] 



(2) 



where An (7) is the largest eigenvalue of the modified 
transfer matrix T . If so inclined, the reader may refer to 
Appendix A to examine the full solution to An (7), orig- 
inally obtained by Derrida and Lebowitz in [29] and an- 
alytically continued by Derrida and Appert in [30] . The 
most important feature of this solution for Aw (7) is that 
it has the scaling behavior 



Aw (7) - ic p 



AT3/2 



(3) 



for fixed filling ratio p/N = p in the scaling limit A*" — > 00 
with "fN 1 / 2 held constant. a p ,b p and c p are all con- 
stants that depend solely on the particle density p. The 
most interesting aspect of this solution is the Derrida- 
Lcbowitz scaling function G, whose form is independent 
of all system parameters. This led [30] to postulate that, 
in fact, G represents a universal scaling function that can 
be used to describe all systems within the KPZ univer- 
sality class. For details about the explicit representation 
of the Derrida-Lcbowitz scaling function, the reader may 
refer to Appendix B. 
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sizes N. More precisely, we suggest that 



FIG. 3: Plot of the left hand side of Eq. (5), rescaled by N 3/2 , 
for systems size N = 2, 4, 6, 8, 10 and 12. Notice how, even in 
the case of very small sizes, these values plotted are almost 
indistinguishable from not only each other, but the adjusted 
DLSF, given in Eq. (7) which they are plotted against. The 
gap in the data occurs due to the slow convergence of the 
plotted solutions (discussed in appendix A) around the point 
B c (as given by Eq. (A3)) 



III. SCALING BEHAVIOR FOR 
INTERMEDIATE N 



Eq. (3) contains a formula for successfully describing 
An (7) in the scaling regime of constant 7./V 1 / 2 in the 
limit of large system sizes. However, when considering 
small systems or large 7 at a fixed system size, Eq. (3) no 
longer adequately describes An (7). Fig. 2 shows just how 
poorly Eq. (3) performs under these conditions. If one 
wanted to extract the scaling function G from numerical 
data for An (7) in order to, e.g., verify its universality in 
a different system, one would have to go to very large 
system sizes N. 

Here, we introduce another approach that can be uti- 
lized to directly verify the scaling behavior of An (7)- To 
this end, we exploit the fact that An (7) not only has a 
well-defined scaling limit as N — ► 00 at constant jN 1 / 2 
(which implies 7^0) but that there is also a well-defined 
scaling limit as N — * 00 at constant 7 for 7 < 0. This 
scaling limit is given [30] by 



A' 



Jim An (7) = Aoo(7) = 



-(1 - cxp[7p])(l - exp[ 7 (l - p))) 
1 - exp[7] 

(4) 

This behavior is not universal, as the explicit dependence 
of the result on the particle density p shows. However, 
the existence of such a scaling limit should still be a uni- 
versal feature of all KPZ systems. 

We propose that a combination of the two scaling limits 
in Eqs. (3) and (4) provides a faithful representation of 
the full function An (7) for 7 < and very small system 



An (7) ~ ^(7) 



apG^N^bp) a^bpf 



N 3 / 2 



24tt 



(5) 



already for very moderately sized N = O(10). This 
agreement is essentially due to the fact that Aoo(7) pro- 
vides an approximation the An (7) behavior that is far 
superior to the linear approximation from which the 
Derrida-Lebowitz scaling function G was originally de- 
rived. 

It is instructive to verify that our new relation Eq. (5) 
indeed simplifies to the two known scaling relations 
Eqs. (3) and (4) in the appropriate limits. For fixed 
7 < the argument of the Derrida-Lebowitz scaling func- 
tion G on the right hand side goes to negative infinity as 
TV — > 00. In this limit it is known [301 that 



3 



24tt 



(6) 



where the sub-leading terms decay exponentially with 
\(3\. Since this leading term in G cancels with the sec- 
ond term on the right hand side of Eq. (5), the right 
hand side vanishes as N — * 00 and Eq. (4) emerges. If 
on the other hand ^N 1 / 2 is held constant as N — > 00, 
7 vanishes. Thus, the function A QO (7) can be expanded 
for small 7. Since Aoo(7) is an odd function, this ex- 
pansion contains only terms with odd powers of 7. The 
linear term precisely yields the term 7c p in Eq. (3). The 
7 term equals the second term on the right hand side 
of Eq. (5) (note that this implies a non-trivial connec- 
tion between the non-universal constants a p , b p , and c p 
and the non- universal function Aoo(7)). All higher order 
terms in 7 become sub-leading as N — > 00 in this scaling 
limit. Thus, Eq. (3) reemerges. As it becomes clear from 
this discussion, the additional term on the right hand side 
of Eq. (11) can either be understood as a 7 s correction 
to the left hand side or as a correction to the DLSF G 
on the right hand side. It is the only term the scaling 
behavior of which allows it to be interpreted as both a 
part of the universal scaling function as well as a part of 
the infinite size solution Aoo(7). For the purpose of this 
paper we will integrate this term into our scaling function 
and call the function 



GW) = Gift + ^ 



(7) 



the adjusted Derrida-Lebowitz scaling function. 

Note that while Eq. (5) has only been shown to hold 
in the two limiting cases discussed above, in actuality, 
it works exceedingly well even for small systems. Fig. 3 
shows the left hand side of Eq. (5), rescaled by N 3 / 2 , 
plotted for N = 2, 4, 6, 8, 10 and 12. The right hand side 
of Eq. (5) is plotted as well, with the appropriate rescal- 
ing factors given by [30]. The quality with which Eq. (5) 
captures the behavior of this intermediate scaling region 
is suprising, especially when considering the small system 
sizes. This empirical observation forms the foundation of 
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FIG. 4: Diagram of our discrete time asymmetric exclusion 
process. The diagram shows an odd time where only particles 
on even sites are considered for hopping. Notice that the 
particle at site 1 cannot move because site 2 is already at 
maximum occupancy. However, for sites 3 and 5, the shaded 
particles can hop to the right with a probability rj, as shown 
by the boxes with dotted outlines. 

the methodology presented herein for numerically under- 
standing the behavior of quantities, such as the particle 
hopping, in systems that currently remain beyond the 
reach of purely analytic methods. 



IV. DISCRETE TIME ASYMMETRIC 
EXCLUSION PROCESS 

In the following, as a specific application of our 
method, we verify the universality of the DLSF for the 
discrete time totally ASEP using the sublattice-parallcl 
updating scheme [13]. In our totally ASEP, we consider 
a periodically bound system of size N where p particles 
can only move to the right, and each site can accommo- 
date up to n particles. During each odd time interval t, 
the positions with odd numbers are evaluated for tran- 
sitions. Particles can hop only if available space exists 
to the right for the particle to move, i.e., only if the site 
does not already contain the maximum number of al- 
lowed particles n. For allowed transitions, particles hop 
with probability rj and stay put with probability 1 — rj. If 
there is more than one particle on a site only one parti- 
cle is considered for hopping. For even time intervals, the 
exact same dynamic occurs at the even numbered posi- 
tions. Notice that, in order to use this sublattice-parallcl 
update scheme, the number of sites N must be even. For 
this system, we solve via a transfer matrix method de- 
scribed in Appendix C for the exponential term, A at, of 
the generating function, Z, of the total particle flux per 
site, yt, in the large time limit, t — > oo. Again, let us first 
define this generating function as 

Z t , N {r,0) = (exp[yy t ]) (8) 

where the brackets (■) denote the ensemble average and 
£1 summarizes the specific parameters contributing to the 




FIG. 5: The appropriately rescaled and adjusted DLSF (in ac- 
cordance with the right hand side of Eq. (13)) plotted against 
the direct measurement of the left hand side of Eq. (13) for 
n = 1, p = p/{nN) = 1/2 and rj = 3/4. The left hand side 
of Eq. (11) has been plotted for N = 6, 10, 14, 18, and 22. 
It becomes difficult to distinguish the values for the larger 
widths from the adjusted DLSF solved in [30] since they lie in 
near perfect agreement. This indicates that the proposed uni- 
versality of the DLSF does hold for this system. The scaling 
factors an and bn in Eq. (11), neither of which depends on the 
finite size effects, are chosen only once for the largest system, 
N — 22, and control the scaling of the adjusted DLSF. 

evaluation of the particle flux per site yt, i.e., 77, p, and 
n. As in the continuous time case examined by Derrida 
and Lebowitz [29], the generating function behaves like 

z t ,N{r,ty ^cx P [\ N { Tl n)t} (9) 

for large times, where A at (7; fi) is the largest eigenvalue of 
the characteristic matrix, Tjv(7;fi), the technical details 
of whose construction are provided in Appendix C. 

Since we will only be able to numerically evaluate 
AAr(7;f2) for relatively small N, we apply the method 
presented in the previous section. Toward this end, we 
introduce the function Aoo(7;f2) describing the infinite 
size behavior of the Ajv(7; f2) as 

A 00 ( 7 ;fi) = lim A A r( 7 ;fi). (10) 

N— *oo 

The new scaling form then becomes 

Aw(7;^) - Aoo(7;^) = j^f 2 ■ (ii) 

For the specific case in which the allowed number of 
particles per site equals one, n = 1, and the system is 
half filled, p = N/2, there exists an exact analytical so- 
lution [16, 17] for A oc (7;f7). 
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FIG. 6: Plot of the right hand side of Eq. (14) versus the 
appropriately fitted left hand side for r\ — 3/4 for various 
N. This plots shows the excellent agreement between the 
calculated data and the Derrida et al. proposed form of the 
scaling function under different parameter choices. The total 
number of particles has been chosen such that p = Nn/2, i.e., 
a half-filled system. For n = 1, N = 22 was used, for n — 2, 
N = 14, for n = 4, TV = 10, and for n = 6, TV = 8. 



This, combined with our numerical transfer matrix 
method for calculating Xn(j;H), allows us to directly 
measure the left hand side of Eq. (11). Since the ad- 
justed DLSF on the right hand side has already been 
solved, only the scaling coefficients an and on are un- 
known. This allows us to numerically calculate the left 
hand side of Eq. (11) and use the results to fit an and 
bn for one (the largest) system size TV. For convenience, 
here we multiply Eq. (11) by TV 3 / 2 

N 3 / 2 (X N (r, fi) - Aoo(7; A)) = a^G^N 1 '^). (13) 

and plot the left hand side of Eq. (13) against (3 = 7 TV 1 / 2 
in Fig. 5 for TV = 6, 10, 14, 18, and 22. The agreement be- 
tween the various curves in Fig. 5 gives a good indication 
that the asymptotic solution for the DLSF indeed applies 
to our discrete time ASEP. 



V. DISCRETE TIME MULTI-PARTICLE 
ASYMMETRIC EXCLUSION PROCESS 

In order to give a further example of the usability of 
our method, we also apply it to the case where we allow 
more than one particle per site. A lattice point with more 
than one allowed state or highway traffic with more than 
one lane both provide very good examples of why such 
a multiple particle ASEP is important in and of itself. 
For n > 1, an analytic solution for Aoo(7;f2) does not 
exist. Thus, we cannot quite follow the method shown 
in the previous sections. However, we can still test both 
the validity of the DLSF, and calculate the non- universal 



scaling factors through other means. Since we have a 
method for obtaining Xn{i] ty, wc take the difference 

Ajv(t;^) - A w _ 2 (7;^) = 

anG(jN^ 2 b n ) _ anG( 7 (A r -2) 1 / 2 b n ) (14) 
A«72 (JV-2)3/ 2 

in order to eliminate the need for the value of A oc (7; f2). 
Notice that the size independent terms (an(70n) 3 / (24tt)) 
cancel out, leaving us with the original DLSF G. Once 
again, we find the scaling coefficients an and on by fitting 
the right hand side to the left hand side of Eq. (14) for 
a single TV. Fig. 6 shows the right hand side of Eq. (14) 
for the largest TV with the numerically fitted values for 
an and on plotted against the left hand side whose value 
is obtained using the same transfer matrix method de- 
scribed in appendix C. While in this approach curves 
obtained for different TV do not overlap each other, we 
note that we obtain similarly good agreement as in Fig. 6 
for small system sizes TV without refitting the scaling pa- 
rameters an and 6n (data not shown). This once again 
results in excellent agreement between our solution for 
the left hand side and the properly rescaled right hand 
side for various parameter values £1 including different 
values for n, 77, and TV. 

Fig. 5 and Fig. 6 each display excellent agreement with 
established KPZ theory and come without any real ad- 
justment or reformulation of existing theory. They are 
simply the direct result of studying our multi-particle 
discrete time ASEP dynamics modeled by a transfer ma- 
trix method. Once an and bn have been determined, 
their values can be utilized in order to numerically com- 
pute the infinite form Aoo(7;f2) from A at (7; 17) for even 
small system sizes by utilizing Eq. (11) — a fact which is 
particularly important in cases where Aoo(7;f2) remains 
unsolved. Of course, given a form for Aoo(7;0), even 
though it be numerical, allows one to calculate proper- 
ties of the system in the thermodynamic limit, such as 
the particle hopping and large deviation function. 

In summary, the fast convergence in TV allows for rapid 
calculation of an entire function Aoo(7;f2), an otherwise 
difficult quantity to compute for most systems. The ac- 
cessibility of Aoo(7; f2) through studies of smaller systems 
has already been put to use in calculating practical quan- 
tities important in sequence alignment [18]. 

VI. CONCLUSION 

In this letter, we have presented a numerical method 
that allows a direct numerical calculation of the universal 
DLSF and it's non-universal scaling factors. In contrast 
to previous numerical studies, our method does not rely 
on sampling but rather on the exact diagonalization of 
sparse modified transfer matrices for very small finite size 
systems. Thus, it is able to capture the full information 
on rare events that the DLSF encodes. As an applica- 
tion, we have extended the universality of the Derrida- 
Lcbowitz scaling function to inherently discrete time hop- 
ping processes with multiple particles per site by directly 
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measuring the scaling function itself. We found that the found in [30] , 
numerically determined scaling function converge toward 

the universal scaling function already for relatively small 2B 1 i/ P (1 — 2B 1 ^ P — \/l — 4B 1 / p ) 

finite systems. Mt) = " " _ ^ _ ^ = 

The method we have outlined is applicable to a variety 

of discrete KPZ systems, including the important open 2B 1 /p 0- ~ 2£? 1 / p + s/l — AB 1 /?) 

boundary ASEP and the open boundary partially ASEP j _|_ i q _ 2B 1 /p + Vl — 4B 1 /?) 
used for modeling true non-equilibrium driven lattice gas 

dynamics. In addition to this, the more complex inter- _ p^ 00 59 (-/Vg — 2)! 

actions of multiple number and types of particles can be 9-1 (piV-iNq — VI ~ 1)! 
modeled using our method. 



(A6) 



7 = f (Ml + 2bW(1 - 2B 1/P - VI - 4P7p)) 
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APPENDIX A: \ N FOR A CONTINUOUS TIME 
TOTALLY ASYMMETRIC EXCLUSION 
PROCESS 



The solution for An (7) can be given by the parametric 
equation [29] 



M7) = -P£B* <^;f_ n| (Al) 



^ (pqy.(Nq-pq-l) 



N J2 B 



(Nq-l)l 
' {pq)\{Nq-pq)\ 



(A2) 



APPENDIX B: THE DERRIDA-LEBOWITZ 
SCALING FUNCTION 

The form of the Derrida-Lebowitz scaling function can 
bc given as 



G{{3) 



,3/2 



Ce~ e de 



30? Jo 1 + Ce- 



2_ [°° 1/a Ce^de 



(Bl) 
(B2) 



for C > — 1. For the region less than j3- = limc^-i /?, 
the analytic continuation of G is given by 



for — B c < B < B c , where the radius of convergence B c 
is given by 



pP(N-pY N -P~> 



(A3) 



This set of equations represent the solution for An (7) in 
the region 7_ < 7 < 7+, where 



^ c! (pq)\(Nq-pq)\ 



9=1 



7+ = -NJ2(-B C 



q=l 



(Nq-l)l 
(pqy.(Nq-pqy. 



(A4) 



(A5) 



Derrida and Appert term this the 'scaling region', i.e., 
the region where Eq. (3) holds true in the limit N — > 00. 

In the strongly negative region, 7 < 7_, the solution is 
given by the analytic continuation of Eqs. (Al) and (A2), 



G{p) = -V^[-M-C)] 3/2 -X>c)v 5/2 ( B3 ) 



q=l 



f 3 = -4^F[-iu(-c)] i/2 -Y^{-c) q q- z/2 m 

0=1 



for < C < 1. 



APPENDIX C: BRIEF DESCRIPTION OF 
MATRIX METHOD 

We shortly describe how we actually calculate 
An (7; fi)- In our brand of the ASEP, we use a sublattice- 
parallel updating scheme where hopping probabilities in 
even and odd time intervals are evaluated separately. Be- 
cause our discrete time processes may be thought of as a 
combination of a number of very simple process occurring 
in some sequential order, we may first examine the base 
dynamic and expand this into the larger picture. Before 
examining the movement at all positions, we study the 
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dynamics of a single hopping transition. We first create 
the transfer matrix T(7 = 0; fi) describing the transition 
probabilities for one pair of sites [17] using particle occu- 
pancy number d as our basis. Next, we modify this ma- 
trix by multiplying all off diagonal elements by the factor 
exp[— 7/iV]. This effectively tags the average number of 
hops per site. For n — 2 particles per site this results in 
the matrix 
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in the basis (00, 01, 02, 10, 11, 12, 20, 21, 22) where 
z = (1 — rf) exp[— 7/iV]. 

From this smaller transfer matrix, we may build up 
this single pair description into the larger iV-site picture 
by taking the tensor product 

N/2 

f N ( r ,n) = <g)T(r,n) (C2) 

k=l 

and eliminating any states that do not contain the right 
number of particles p. This gives us the matrix that 



models particle hopping from sites in one time interval. 
Converting this to the appropriate basis for the next time 
interval can be done by utilizing the translation operator 
C defined such that 

C\d di . . . d N -i) = \did 2 . . . djv-ido)- (C3) 

After exploiting translational invariance and up-down 
mirror symmetry in order to reduce the size of our 
state space we obtain the additional identity C = C _1 
on this reduced state space. Then the matrix product 
f even f odd = f N (C- x f N C) = (f N C) 2 describes the par- 
ticle hopping of our discrete time ASEP for all N sites. 
These dynamics can be viewed as a Markov process on a 
(11 + 1 ^-dimensional state space of the equal time dif- 
ference vector \d(0,t),d(l,t), . . . ,d(N,t)). Solving for the 
largest eigenvalue £w(7; of TWC gives us the particle 
hopping function for finite size [16, 17] 

\ N ( T ,n) = ln6r( 7 ;fi). (C4) 

This is the function used in order to produce the re- 
sults given above. Calculating this largest eigenvalue 
is somewhat challenging, since this matrix description 
grows very quickly with TV. However, the matrices are 
very sparse with the number of non-zero matrix elements 
growing almost linearly with the matrix dimension. This 
makes it possible to numerically obtain the largest eigen- 
value required in Eq. (C4) using the implicitly restarted 
Arnoldi method [34] for matrix dimensions up to around 
10 5 . This allows us to plot the entire function Xn(^;H,) 
in very reasonable timescales. 
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